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Abstract. As a stochastic model for quantum mechanics we present a stationary 
quantum Markov process for the time evolution of the Wigner function on a lattice 
phase space Zj\j x Zj\j with N odd. By introducing a phase factor extension to the phase 
space, each particle can be treated independently. This is an improvement on earlier 
methods that require the whole distribution function to determine the evolution of a 
constituent particle. The process has branching and vanishing points, though a finite 
time interval can be maintained between the branchings. The procedure to perform a 
simulation using the process is presented. 
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1. Introduction 

Various studies have been made of quantum mechanics by using trajectories. For 
example, in Feynman's path integral method, the transition amplitude between the 
initial and final points is calculated by adding all the phase factors determined from 
the trajectories connecting these points [1]. In this case, the trajectories are interfering 
alternatives and not those of observed particles. The method proposed by Bohm is 
another example. In this method, the trajectories are treated as if they were observed 
and are calculated deterministically [2] . These methods give alternative formulations of 
quantum mechanics. 

The stochastic approach is another method. It is used not only for the formulation of 
quantum mechanics, but also for practical calculations. The use of a stochastic method 
for practical calculations is generally called the quantum Monte Carlo method (QMC) 
and is widely used in the simulation of quantum systems. Typical algorithms in the 
QMC are the random walk Monte Carlo method (RWMC), the diffusion Monte Carlo 
method (DMC), and the path integral Monte Carlo method (PIMC). These methods are 
efficient ways of determining the wave function and energy of a ground state [3,4]. They 
are based on the formal analogy that the Schrodingcr equation becomes the diffusion 
equation when t is substituted by an imaginary time ir. 

A solution of the diffusion equation can be constructed using Brownian motion. 
We can generate sample paths of Brownian particles by a random number generator 
and can follow the paths independently as Brownian motion is a stationary Markov 
process, which means that the diffusion equation can be solved by a stochastic method 
using a random number generator. Many studies have been conducted on the diffusion 
equation and Brownian motion. It is known that a stochastic estimation is generally 
more efficient than direct calculations when the system has many degrees of freedom. 
In direct methods, the amount of calculation required generally increases exponentially 
with the number of degrees of freedom D [5] . On the other hand, for stochastic methods 
the amount of calculation is of order D for one trial, though the estimated error is of 
order 1 / ^/n where n is the number of trials. Very naively, the Schrodinger equation is 
a diffusion equation along an imaginary time axis. However, the quantum stationary 
Markov process for the Schrodinger equation, which corresponds to Brownian motion for 
the diffusion equation, has not yet been established. Formulating an analogous process 
to Brownian motion in quantum mechanics would be useful for estimating the time 
development of wave functions, especially in higher degrees of freedom. 

Because of the probabilistic nature of quantum mechanics, several stochastic 
formulations of quantum mechanics have been investigated [6,7]. A typical example 
of such a stochastic process was formulated by Nelson in 1966 [6]. In this method, the 
behavior of the constituent particles is determined by a stochastic differential equation 
and their distribution function coincides with the absolute square of the wave function. 
In this sense, the trajectories in this method are those of the observed particles. We 
can apply this method to estimate the tunneling time by making simulations of the 
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trajectories [8]. However, the drift term in the stochastic equation is calculated from 
the wave function, and hence the Schrodingcr equation must be solved to construct 
the process. A general method to construct a stochastic process without solving the 
Schrodinger equation was developed by Guerra and Marra [9]. In their method, the 
total distribution of particles is required to determine the next step of a constituent 
particle, though obtained process is a Markov process. It is possible to determine the 
time evolution of all particles simultaneously, but each individual particle cannot be 
treated independently. These methods treat the distribution function directly, i.e. the 
absolute square of a wave function. 

A way to avoid this difficulty is to search for a stochastic process for the time 
evolution equations which are equivalent to the Schrodinger equation. One such equation 
is the quantum Liouville equation for the Wigner function. In this paper, we discuss a 
stationary quantum stochastic process for the Wigner function in a lattice phase space. 
A similar approach was considered by Cohendet et al. [10]. They constructed the 
Wigner function on a lattice constituting of an odd number of lattice points, and found 
a background quantum Markov process. 

The Wigner function is defined by Wigner as 



in phase space, where is a wave function [11,12]. The Wigner function gives the 
density distribution in the configuration space if it is integrated along the momentum 
direction: 



and gives the density distribution in the momentum space if integrated along the 
configuration direction: 



This property is the so-called marginality of the Wigner function. In this sense, the 
Wigner function can be regarded as a kind of distribution function in phase space and 
plays an important role in quantum mechanics, especially in the quantum-classical 
correspondence. However, the value of the Wigner function can be negative, so that it 
is not a distribution function itself. 

Cohendet et al. regarded a lattice with odd N lattice points as a cychc group Zjv of 
order N and constructed a Wigner function in the lattice phase space Zjv x Z^, i.e. the 
direct product of Zj^ and its dual (= Zj^). The Wigner function can also be negative on 
a finite lattice, but the value is bounded. They extended the phase space Zjq x Zjsi by a 
dichotomic variable a of value d=l and considered a new function on the extended space 
adding a constant to the Wigner function. The function is normalized and positive 
as the original function is bounded, so it can be regarded as a distribution function of 
particles in the new space. They derived the time evolution equation for the distribution 
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function from the quantum Liouville equation, which determines the time evolution of 
the Wigner function, and constructed the background stochastic process of the particles 
using the Guerra-Marra method. Although the obtained process has a Markov property, 
the whole distribution function is required to determine the evolution of a constituent 
particle. The transition probability is time dependent and not stationary. 

In this paper we introduce a stationary Markov process that enables us to treat all 
particles independently, by extending the phase space by a phase factor U{1) on an odd 
lattice. The process has branching and vanishing points and a finite time interval can 
be maintained between the branchings. 

2. Stationciry quantum Mcirkov process for the Wigner function 

2.1. Wigner function on a lattice 

Let be an odd number. We define N x N phase and shift matrices Q,P as 



Q — {Qnm) = {^"'^n,m) — 



UJ 



V 

/ 
1 



_iv I 1 
UJ 2^2 



(4) 



P — {Pnm) = (^n+l,m) 



1 \ 



1 

1 



V 



(5) 



1 0/ 



where m.n G Zj^f and are labeled from (A^ — l)/2 to — (A^ — l)/2 with step -1, d is the 
Kronecker delta function modulo A^, a; is a primitive A^-th root of unity, e.g. exp{27ri/N). 
The matrices Q, P correspond to the phase and shift operators 



(6) 



respectively, in continuous space-time, where a is a lattice constant. The Weyl matrices 
W{m, n) and Fano matrices A(m, n) are defined as 



-2m 



(7) 



A(m, n) = W{m,n)T, 



(8) 
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where T is a matrix given by 

/ 1 \ 



V 1 / 

The Weyl matrices W{m,n) can also be written as 

W{m, n) = uj'^rnnp-2mQ2n ^ p-mQ2np-m ^ Qnp-2mQn 



(9) 



(10) 



The A^^ Weyl matrices and Fano matrices are both complete and orthonormal in 
the trace norm: 

Tr[W\m,n)W{m,n)] = N5m,n. 

T:[At(m,n)A(m,n)] = N5m,n, 

W\m, n) = W{—m, —n), 

A"^(m, n) = A(m, n), (11) 

so that arbitrary matrices can be expressed as a linear combination of these matrices. 
We call the coefficients that appear in the Weyl and Fano expansions the Weyl and Fano 
coefficients, respectively. 

In the continuous case, the Weyl operator is defined as 



where 



Qc 



Pc = e'". 



The Wigner function is written as 

Pciq,p) = ^Tr[pW,(q,p)T,] 



(12) 



(13) 



(14) 



where the operator Tc is given by 

/oo 
dq\-q){q\. (15) 
-oo 

The Wigner function on a lattice is defined in an analogous way to the continuous case: 

(16) 

It is easily seen that this definition satisfies the marginality condition on the lattice. 



p(m,n) = ^Tr[pW{m,n)T] = lT:[pA(m, n)]. 
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2.2. Quantum Liouville equation for the Wigner function 
Let H{m, n) be the Weyl coefficients of the Hamiltonian, 

H ^Y.^('^^^)Hijn,n), (17) 

m,n 

and p(m, n) be the Fano coefficients of the density matrix, i.e. the Winger function, 

p = ^ A(m, n)p(m, n). (18) 

m,n 

The quantum Liouville equation 

ihj^p = [H, p] (19) 
can be written using these coefficients as 

-^p(m, n) = ^Y. {H{m -m',n- ^')a;-2("»"'-""^') 

CLv Civ / / 

m' 

- H*{m -m',n- n')a;^^"'"'~""''^}p(m', n'). (20) 

Wc modify the equation as follows for later convenience. We perform a polar 
decomposition of the Weyl coefficients of the Hamiltonian: 

H(m, n) = h(m, n)a;^("^'"), (21) 

where ^(0, 0) is from the Hermiticity of the Hamiltonian. The time evolution equation 
for the Wigner function (20) can be rewritten by this decomposition as follows: 

— p(m, n) — — ^2 h{m — m', n — n') 

^ m',n' 

^ ^^-2(mn'-nm')+e(m-m\n-n')-N/4: ^2{mn'-nm')-e{m-m',n-n')+N/4^^^^/ ^/-j ^22^ 

The term in the right-hand side equals zero when m — m',n — n', so that it can be 
omitted in the sum. 

2.3. Generator of the Markov process and time evolution equation for (p 

We extend the phase space Zn x with phase factor U{1) and consider a non negative 
function (p{a, m, n) on C/(l) x Zj^ x Zj^ that has the following relation with the Wigner 
function: 

p{m,n) = / dae^°'ip{a,m,n). (23) 



Since the Wigner function is bounded in modulus by 1/A^ on a lattice 
1,1 
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the function (p{a, m, n) exists, although it is not uniquely determined. We can take 

(f{a,m,n) = Lp{—a,m,n) (24) 

without loss of generality as p{m, n) is real. 

We consider the following time evolution equation for (p{a,m,n) within a short 
time period At: 

(p{a,m,n;t + At) — J da' M/^t{a,m,n;a',m',n')(p{a',m',n';t) (25) 

m' ,n' 

where MAticx, m, n; a', m', n') is given by 

MAt{a, m, n; a', m', n') = (1 - DAt)S{a - fAt{a'))Sm,m'^n,n' 

H —h[m — m ,n — n ) 

x{S{a -a' - 27r(-2(mn' - nm') + e{m -m',n- n') - N/A)/N) 

+ 5{a -a' - 27r(+2(mn' - nm') - e{m -m',n- n') + N/A) /N)}, (26) 



with 



\h\ = Y.h{m,n). (27) 



We assume the function fAt{(^) in this equation satisfies 

(1 - DAt) cos /a* (a) = cos a, (28) 

^mjAtia)=a. (29) 

The definition of MAt{<y, "m-, n; a', m', n') shows that particles with ratio DAt move 
to other phase points with probability h{m — m',n — n')/\h\ within At. On the other 
hand, other particles stay at the same phase points and change their phase from a 
to /At(tt) deterministically. This means that M^t represents the time evolution of a 
Markov process for (p{a, m, n). 

2.4- Time evolution of the Wigner function 

Wc consider the time evolution of the Wigner function p{m, n) when <^(q;, m, n) evolves 
following Eq. (25): 

p{m,n]t -\- At) — j da j da' ^ e^°'MAt{oi,m,n;a',m',n')(p{a',m',n';t) 

m',n' 

— p(m, n; t) + [ da' V] h{m — m', n — n') 

2|/i| 

X (g^"'^'~2(mr!.'— nm')+6(m— m',n-n')-Ar/4 
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We have 



-^{p(m, n; t + At) - p{m, n; t)} 
— — — ^ h{m — m' , n — n') 

m',n' 



^ ^^-2{mn' -nm')+e{m-m' ,n-n')-N / i _|_ ^+2(mn'-nm')-6'(m-m',n-n')+A''/4^ 

X da'e''''ip(a',m',n';t) (31) 
and in the At ^ hmit 

-^p(m, n) = -2- ^ /i(m - m', n - n') 

^ ^^-2(mn'-nmO+6'(m-m^n-n')-A^/4 _|_ ^+2(mn'-nm')-6'(m-m',n-n')+Ar/4-^^^^/ ^/^ ^^2) 

This shows that if /^{a, m, n) evolves following the time evolution equation (25) and 
we take 

D/2\h\ = 1/h, (33) 

the function p(m, n) obtained from the distribution function (/?(q;, m, n) satisfies the 
quantum Liouville equation for the Wigner function. 

2. 5. Generator of the Markov process for (f 

The generator of a Markov process 

A{a, m, n; a', m', n') is generally defined by 

A{a,m,n\a'm\n') = lim — ^ . (34) 

^ ' ' ' ' ' ' At->o At 

The generator corresponding to MAt(Q;, m, n; a\ m', n') is given by 

A{a, m, n; a', m', n') = ^ ^ 8\a - a')5m,m'5n,n' 

-DS{a - a')Sm,m'8n,n' + ^^("^ - m' , u - n') 

X {5{a -a' - 27r(-2(mn' - nm') + e{m -m',n- n') - N/A) /N) 

+ 5{a -a' - 27r(+2(mn' - nm') - 9{m -m',n- n') + N/A) /N)}. (35) 

This generator is singular at a = 0, tt. These singularities come from the fact that 
when we set 

/At(tt) = 01 + Act (36) 
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in Eq. (28) and Eq. (29), we have 

At/Aa = --^tana = (37) 

and At becomes zero at these points. This means we cannot take the time evolution 

for any choice of Act at these points. To avoid this difficuhy, we have to choose the 
initial ip{a,m.n) to be zero at a = 0,7r and perform branching without changing the 
contribution to the Wigner function when a particle reaches these points. For example, 
a particle which reaches a = can be replaced by two particles at a = +7r/3 and 
a — — 7r/3. 

In the limit Af — > 0, the time evolution equation for ip{a,m,n) satisfies 

. . d D \ / \ 

o?(q;, m, n) — — 09(0;, m, n) — D(p{a, m, n) 

oa tan a 

f D ~ 

+ / da' ^ — —h{m — m', n — n') 

m',n' 2|^| 

x{S{a -a' - 27r(-2(mn' - nm') + e{m -m',n- n') - N/A)/N) 

+5{a -a' - 27r(+2(mn' - nm!) - 0{m -m',n- n') + A^/4) /N)} 

X ip{a',m',n'). (38) 

The quantum Liouville equation for the Wigner function p(m, n) also can be derived 
directly from this equation. We observe that Eq. (38) reduces to Eq. (22) if the part 

d D 



da tan a 



(p{a,m,n) — Dip(a,m,n) (39) 



vanishes after the integration with respect to a multiplying e*". Paying attention to 
(fi{a) — (fi{—a), the integration 



2^ . d 1 

e'"- (p(a)da 

oa tan a 



becomes 



r d I , , , . / ^ rdcosa i , , , 

! / cosctT" (piajda = 2\sma(p{a)L — 2 / — (piajda 

Jo ao; tan a Jo oa tana 

— 2 cosaip{a)da = / e'^'^Lp{a)da (40) 
Jo Jo 



and Eq. (22) holds. 



3. Construction procedure 

In the following we summarize the procedure to determine the time evolution of 
the Wigner function from the Markov process of particles with distribution function 
(fi{a, m, n). 
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3.1. Preparation 

1-1) Prepare an additional space U{1) on each phase space point {m,n). 



1-2) Calculate the Wigner function from the initial wave function. 

1-3) Expand the Wigner function into the distribution function and make an appropri- 
ate choice of (p{a, m, n). 

1-4) Perform the Weyl expansion of the Hamiltonian and calculate the Weyl coefficients 
H{m, n). 

1- 5) Perform the polar decomposition of H{m,n) and calculate 9{m,n). 

3.2. Rules of the process 

2- 1) Particles of ratio DAi move to other phase points within At. Others stay at the 
same phase point. 

2-2) The probability P(m, n\ m', n') for moving from (m', n') to (m, n) is determined by 
P(m, n; m', n') — h{m — m', n — n')/\h\, (41) 
with Eq. (27). 

2- 3) The phase change is deterministic for particles without jumps to other phase points. 
The phase is rotated from a to /At(Q;). 

3.3. Simulation 

3- 1) Prepare an ensemble of particles with the initial distribution function. 

3-2) Make transitions for each particle independently in the extended phase space fol- 
lowing the probability rule above within a period of time. 

3-3) For a jump to another phase point, add the phases ;^-2{'^n'-nm')+6'(m-m',n-n')-Af/4 
and uj'^irnn'-nm')-e{m-m.' ,n-n')+N/A ^-^j^ probabifity 1/2. 

3-4) Pull back the particle when it reaches the stopping point a = 0, tt with branching. 

3-5) After performing the process for each particle, integrate the phase part multiplying 
e*" to the density distribution function (^(a, m, n). The resultant function is the Wigner 
function p{m,n) under Eq. (33). 
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The ensemble must be remade if the number of particles increases. 

Wc show an example foT N — 3 just as a clue to applications. When the Hamiltonian 
and initial wave function are given by 



H = S, 



/ 1 \ 





( 1 \ 

V 1 / 



with 



where 



< >= 1, < Sy > = < Sz >= 0, 



/ 1 \ 



10 1 
1 



/ -1 \ 



- — 



1 -1 
1 



the Wigner function and an extended distribution function are given by 



p(m, n) 



1 - V2 

12 



, ip{a,m,n) 



V2-1 



12 



(S{a-—) + 6{a + —)j, 



for (m,n)=(l,l),(l,-l),(-l,l),(-l,-l), 

l + 2x/2 



p{m, n) 



12 



for (m,n)=(0,l),(0,-l), 

p(m,n) = 1/12, ip{a,m,n) 
for (m,n)=(l,0),(-l,0). 



1 + 2^ 



12 



1 / TT \ 

p(m, n) = 1/3, (/^(q;, m, n) = \^6{a - g ) + (^(a + g ) j , 

for (m, n)=(0,0). The Weyl coefficents have only two non-zero values, i.e., 

^(0, 1) = i/V3, HiO, -1) = -i/V3, 
hence the stochastic process is simply determined from 

P(0,l) = P(0,-l) = l/2, 

^(0,1) = ^(0,-1) = 3/2. 



4. Summary 

If we consider a lattice constituting of N odd lattice points as a cyclic group of 
order N, the corresponding phase space is x Z^. By extending the phase space 
by a phase factor f/(l), we can construct a stationary Markov process in the space 
U{1) X Zn X Zn- This process has the property that multiplying the phase factor to 
the density distribution function of the process and integrating the phase part gives a 
function which satisfies the quantum Liouville equation for the Wigner function. 
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There exist branching points at the phase 0, tt, but we can maintain a finite time 
interval before the branching by making an appropriate choice of the initial condition 
and branching method. There is no contribution from the distribution at the phase 
±7r/2. The process does not need to be continued for particles reaching these points 
as the stochastic rule is symmetric around these points and they give no contribution 
at future times. In this sense, these points can be treated as vanishing points. The 
estimation of branching and vanishing ratio is a problem for future studies. 
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